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Abstract 

We present a mathematica package that performs the symbohc calculation 
of integrals of the form 



where and jfj,{x) denote spherical Bessel functions of integer orders, 

with u > and /jl > 0. With the real parameter u > and the integer 
n, convergence of the integral requires that n + u + /j, > 0. The package 
provides analytical result for the integral in its most simplified form. The 
novel symbolic method employed enables the calculation of a large number of 
integrals of the above form in a fraction of the time required for conventional 
numerical and Mathematica based brute-force methods. We test the accuracy 
of such analytical expressions by comparing the results with their numerical 
counterparts. 
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Licensing provisions: 

Programming language: Mathematica 7.1 

Computer: Any computer running Mathematica 6.0 and later versions 
Operating system: Windows Xp, Linux/ Unix 
RAM: 256 Mb 
Classification: 5 

Nature of problem: Integration, both analytical and numerical, of products of two 
spherical bessel functions with an exponential and polynomial multiplying factor 
can be a very complex task depending on the orders of the spherical bessel func- 
tions. The Mathematica package discussed in this paper solves this problem using 
a novel symbolic approach. 

Solution method: The problem is first cast into a related limit problem which can 
be broken into two related subproblems involving exponential and exponential in- 
tegral functions. Solving the cores of each subproblem symbolically sets the stage 
for an involved expression tree parsing and manipulation to obtain the most sim- 
plified analytic expression for the initial problem. 

Running time: 1 min for typical values of the arguments and can be several mins 
for large values of the input variables. For the test data included, about an hour. 

PACS: 02.30.GP; 02.60.Jh 

Keywords: Symbolic integration; Spherical Bessel function 
1. Introduction 

Integrals of products of spherical bessel functions multiplied by an ex- 
ponential times a polynomial occur in a variety of fields. Nuclear physics, 
scattering theory, calculation of Feynman diagrams, hydrodynamics and elas- 
ticity theory are notable examples where such integrals appear [H EJ IH EJ 
El [71 m [9]. Several authors have proposed various numerical modules that 
can be used to obtain accurate numerical integrations of these functions 
pn [TU] . While most of these deal with general Bessel functions, spherical 
bessel functions actually appear in most physical problems that have spher- 
ical symmetry. A case in point is scattering problems. 

Recently, the application of the density matrix expansion [121 [13] to the 
HF energy obtained from chiral EFT two- and three-nucleon interactions at 
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N^LO \T5\ [T7t [T5] required analytical expressions for integrals of the form 

POO 

I{u,n,fi,u) = / e~''^'^x"'ju{x)j^{x)dx (2) 

where ju{x) and j^X^) denote spherical Bessel functions of integer orders, 
with u > and yU > 0. The real parameter u is such that u > and n 
is an integer constant. Using the asymptotic properties of spherical bessel 
functions and convergence requirement as the argument x goes to zero, one 
obtains the constraint n + u + fi > 0. 

In addition to needing analytical expressions of I{u, n, fi, z/), the applica- 
tion of Ref. [H] required the evaluation of several hundred of such integrals. 
Furthermore, in a large number of cases, the direct implementation of these 
integrals in Mathematica fails to give closed expressions. Finally, it should be 
mentioned that as the order of the spherical bessel functions u and fi as well 
as the exponent n grow, the expressions one obtains for the integrals from 
the direct implementation become unwieldy if and when one has a closed ex- 
pression for those integrals. All these problems necessitated the development 
of a dedicated Mathematica package that is the subject of the present paper. 
The paper is structured as follows: in the next section we discuss the main 
symbolic algorithm and its limitations. This is followed by section |3] that 
presents its Mathematica implementation. Section |4] is devoted to actual 
details of the package and to an illustrative example that shows how one can 
use the developed package. The last two sections contain the conclusions and 
the appendix. 

2. Symbolic Integration 

In this paper, the term symbolic integration signifies a method of com- 
puting complicated integrals using the techniques of expression tree parsing, 
extraction, deletion and replacement with the aim of obtaining a simplified 
expression. The final expression must be the same as the one obtained ana- 
lytically, if possible, whereas for given input values, it should evaluate to the 
same value as the corresponding numerical integration. 

2.1. Recasting the problem 

The starting point for the symbolic integration algorithm is the definition 

rv 

F{y,u,n,fi,u) = / e-''/''x'^j^{x)j^,{x)dx , (3) 
Jo 
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that allows the rewriting of Eq. ([2]) as 

I^UyriyU, fj.) = Urn F{y,u,n, fi,^) , (4) 

y~*oo 

which is valid for all convergent integrals. The essence of such a simple step 
relates to the fact that Mathematica can evaluate F{y, u, n, /i, u) for all cases 
of interest (note the specified domain). One can separate the resulting ex- 
pression for F{y, u, n, fi, v) into a part that contains an explicit y dependence 
and another part without y dependence, viz, 

F{y, u, n, /i, z/) = G{y, u, n, /i, u) + H{u, n, fi, u) . (5) 

Since the most troublesome part is the evaluation of the limit as y goes 
to infinity, the next several steps pertain to the symbolic evaluation of this 
limit. Direct evaluation of the limit of G{y, u, n, fi, v) as y goes to infinitydoes 
not seem to work or takes an impractical amount of time for a large part of 
the parameter space (n, /x, v). This is what necessitated the development of 
the core of the symbolic algorithm. 

2.2. The core of the algorithm 

Extensive tests performed for a large number of combinations of parame- 
ters n, /i, V show that Mathematica always expresses G{y, u, n, fi, v) in terms 
of the exponential integral function and exponential functions of the form 
Qf{y,u) -y^^j-^ere f{y,u) is some function. The exponential integral function is 
defined formally as [2l I16j^ 



Ei{x) = - 




Furthermore, one finds that the y-dependence of G{y, u, n, /i, u) is completely 
absorbed by f{y,u) and by the exponential integral functions. Additionally, 
G{y, u, n, /i, u) does not contain products of exponential functions and expo- 
nential integral functions. Hence, the generic expression that Mathematica 
produces for G{y, u, n, /x, z/) can be written as 

N M 

G{y, u, n, = Y, a*(M)e^'(^'") + ^ bj{u)Et{gj{y, u)) , (7) 
i=i j=i 



^ There is a sign difference between the definitions given in Rcf. [2 and Ref. [16]. That 
of Ref. [16] is adopted in this work. 



4 



where and M are integers whereas aj('u), fi{y,u), bj{u), gj{y,u) denote 
some unspecified auxihary functions. At this point, the algorithm does not 
require the exact form of these auxihary functions to be specified. FinaUy, we 
note that the structure of G{y, u, n, /x, u) as given by Eq. ([T]) can be understood 
adequately by applying the method of partial integration to Eq.(|3]). 

In the evaluation of the limit of G{y, u, n, fi, u) as y goes to infinity, the 
first group of terms terms in Eq.Q, viz, the ones that contain exponential 
functions can be simplified directly. Thus the problem narrows down to 
the evaluation of the limit of the second set of terms in Eq.Q, specifically 
the ones containing exponential integral functions. In line with that, the 
following observations are integral parts of the algorithm: 

(i) There appear only three distinct exponential integral functions, i.e. 
M = 3. Their argument functions, gj{y,u), have the form 



giiy,u) 

92{y,u) 

93{y,u) 



u 

-i2y 



y 

u 



i2y 



u 



(8) 



(ii) The real parts of the hj{u) coefficient functions satisfy the constraints 



Re Y^h, 

Re[h2{u)\ = Relhiu 
(iii) The imaginary parts of the bj{u) coefficient functions satisfy 

= 0, 



(9) 



Im [bi{u)~\ 



-Im [63 ("U 



(10) 



where Re and Im extract the real and imaginary parts of the arguments 
respectively. As discussed in the implementation section (see section [3]), all 
these conjectures can be verified during the program's progress. Hence, in 
the evaluation of the limit of G{y, u, n, /j,, v) as y goes to infinity and after all 
the symbolic manipulations that need to be performed are marked out, one 
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is left with the evaluation of the limits 



hi{u) 
h2{u) 



lim 
lim 



Eil-i2y 



Ei i2y 



- Eili2y 



2Ei 



y_ 

u 



, (11) 
(12) 



Analytic evaluation of these limits as outlined in the appendix 

and h2{u) = —i2TT. 



tion 



7.1) yield the results hi{u) 



see see- 
As mentioned 



in section |3| Mathematica expectedly gives the same answer though, in a 
very complicated form as demonstrated in Eqs. (14) and (15). Refer to sec- 



tion [3^ for the role of the simplifier part of the Mathematica package in this 
context. The subsequent steps of the algorithm comprise making use of the 
listed conjectures and a host of symbolic manipulation operations to obtain 
the most simplified expression for I{u, n, fi, u). 

2.3. Comments on the algorithm 

The symbolic integration algorithm given in the previous section is mostly 
based on Mathematica "symbolic experiments" using a large number of val- 
ues for the input variables n, /i and v. As such, most part of the algorithm 
is based on empirical observation. Even so, all of the conjectures in the algo- 
rithm are verifiable at each stage of the symbolic computation, as discussed 
in section |3j It would of course be beneficial to verify those conjectures ana- 
lytically as our intensive and systematic testing on a vast number of values 
for the input variables did not locate a single exception to these conjectures. 



3. Mathematica Implementation 

The Mathematica implementation of the symbolic algorithm consists of 
three main parts: parser, simplifier and checker. In the following, we discuss 
the role and the implementation of each part. 

3.1. Expression parser 

The first task of the expression parser is to take u, n, fi, u) as an 
input and extract G{y,u,n, fi.u) and H{y,u,n, ^.u). Hence, we define the 
input E{y,u,n, ^.u) which we call Fint as 

Fint = Expand[Integrate[Exp[-x/u] n SphericalBesselJ[fi, X ] 

SphericalBesselJ [v , x], Assumptions-^{lm[u]==0, u > 0}]], 
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where the expression is expanded to make the expression tree available for 
the subsequent parsing steps. The key Mathematica functions to extract 
G{y,u,n, n.u) and H{y,u,n, fi.u) from Pint are Position and Head, along 
with a few other helper functions. 

In the next several steps, all the conjectures mentioned in section |2] are 
verified by parsing and traversing the expression tree. In connection with this, 
the Mathematica function Coefficient turns out to be handy. Specifically, the 
conjectures that need verification are (i) the existence of only exponential and 
exponential integral functions in G{y, u, n, n, v) (ii) that the terms containing 
the exponential function vanishe in the limit of y going to infinity (iii) the 
set of conditions specified in Eqs. (|8]), ^ and (10). The expression parser 



aborts the calculation and produces error message if any one of the above 
conjectures fails to hold. As a side note, we mention that for the domain 
specified in section [T| all the conjectures have been verified (see section 



4.3). In the following step, the expression parser makes use of the two basic 



limits identified in Eqs.([l9|) and (20) to rewrite I{u,n, ^, u) as 



I{u, n, /i, z/) = H{u, n, fi, v) + Ci(u, n, /x, z/)/ii(n) + C2(m, n, /x, h')h2{u) , (13) 

where the coefficient functions Ci(m, n, fi, v) and C2(m, n, /i, z/) are determined 
using the Coefficient Mathematica functions and some additional parsing. 



The corresponding Mathematica code for Eq.(13) reads 



lint = HIint + Coefl*Liml + Coef2*Lim2 , 

where lint is the Mathematica variable for I{u,n,fi, u), HIint is the Mathe- 
matica variable for H{u, n, /i, z/), Coefl is that of Ci{y, u, n, fx, z/), Liml is for 
hi{u), Coef2 is that of C2{y,u,n, fi, v) and Lim2 is for h2{u). 

Paving the way for the discussion of the simplifier part of the package, 
note that Mathematica produces hi{u) and h2{u), prior to any simplification, 
under the form 
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and 



hi{u) = 2(1 + 4^2)2 Log[-22 - l/u] + Su' hog[-2t - 1/u] 

+I6u^ Log[-2i - l/u] + Log[2i - l/u] + 16m^ Log[2i - l/u] 
- Log[-2z + l/u]- Log[-2z + l/u]- IQu^ Log[-2z + 1 /u] 

+8u^ Log[- — — ] - 8u^ Log[2i + l/u]- IQu^ Log[2i + l/u] 

+8m^ Log[-l - 2m] + IQu^ Log[-l - 2m] + 8u^ Log[-l + 2m] 

111 

+ 16m^ Log[-l + 2m] - IQu^ LogM - 32m^ LogM + Log[ ] 

i — 2u 

-Log|2> + Un] + Log|— 1 - Log[— ^1 - Log|-— ] 

4 r 1 , o 2 - r iu , r iu 



+ Log[2z-lHV (14) 



hiu) = 2(1 ^4^2)2 (8^'lQg[-2^ - + 16nMog[-2z - l/u] 

- log[2i - 1/m] - 8^2 log[2i -l/u]- IQu^ \og[2i - 1 /u] 
+ log[-2i + l/u] + 8u^ log[-2i + l/u] + 16M^log[-2i + l/u] 
-log[2i + 1/m] - 8uHog[2i + l/u]- lewMogpi + l/u] 
+8m^ log[-l - 2m] + 16u^ log[-l - 2m] - 8^^ log[-l + 2m] 

tU 111 

-IQu^ log[-l + 2m] + log[, ] + 8^2 log[ 



i — 2u i — 2u 

+16„Mog[— 1 - log[-^l + logl-— 1 

— log[r — — ] — 8u iog[-. — — ] — 16u log 



i + 2u i + 2u i + 2u 

+log[-2z - 1/u]^ , (15) 

which are far more complex than the simple expressions that obtained from 
the ouputs of the simplifier or analytic derivation as discussed in appendix 
EH 
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To summarize, the key task of the the expression parser is to implement 



the Mathematica code outhned just after Eq.(13), thereby arriving at a com- 
pletely analytical but unnecessarily complicated expression for J(m, n, /i, z/). 
Details of the implementation can be found in the submitted Mathematica 
code. 

3.2. Expression simplifier 

At this point, we have at hand analytic expressions for I{u, n, /x, u). How- 
ever, as can be expected from the previous discussion, the expressions are 
extremely complicated and barely useful. To illustrate, the expression one 
obtains for I{u, —1, 1, 1) is 



/(m, -1,1,1) = ' «J , „J , &[ „J 



+ 



4 24m2 I92ti4 16^2 12m 

Log[2^-i] Log[22-i] 2Log[2z-i] Log[-i] 

192^4 16u2 12u 96^4 

Log[-i] ^ Log[-n] , Log[-u] ^og[^^] 



8^2 96^4 8n2 192^4 

Log[3T^] ^Log[^] Log[^] Log[^] 



16n2 12u 192^4 16^2 



^Log[-i+2... Log[^:^] - Log, 

12m ^^i-2u 'i + 2u 

+ 8^2 Log[-2i - -] 
u 

i . 1, „ . 1. 



■ Log[-2z - -] + 8u^ Log[-2? - 



u 

1. _ . 1, o_ . 1 



12m(1 + 4^2)2 

+ IQu^ Log[-2i - -1 - Log[2z - -1 - 8^^ Log[2i 

u u 

- IQu^ Log[2i - -] + Log[-2z + -] + 8^^ Log[-2i + -] 

u u u 

+ 16m^ Log[-2i + -1 - Log[2i + -1 - 8^^ Log[2i + -1 

u u u 

- l&u^ Log[2i + ^] + Log[-l - 2iu] + IQu" Log[-l - 2iu] 



- 8m^ Log[-l + 2iu] - 16m^ Log[-l + 2iu] + Log[ 



i — 2u 
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+ Log[^— ] + 16ii^ Log[-— ] - Log[- 



i — 2u i — 2u —i + 2u 



+ Log[2i - -] + Su^ Log[2i - -] + 16u^ Log[2i - -] 

w 

- Log[-2i + -] - 8u'^ Log[-2i + -] - 16w^ Log[-2i + -1 

u u u 

11 1 

- Log[2i + -] - S-u^ Log[2i + -] - 16m^ Log[2i + -1 

u u u 

+ Log[-l - 2iu] + 16m^ Log[-l - 2iM] + 8u^ Log[-l + 2w] 
+ 16ii^ Log[-l + 2iii] - 16u^ Log[u] - 32k^ Log^ 

+ 8u^ Log[- — — ] + 16^i^ Log[, — — ] - Log[ 



i — 2u i — 2u —i + 2u 



- Log[-— -] + Log[— -] + Log[-. 



i + 2u i + 2u i + 2u 



It should be noted that, the above expression is obtained for relatively 
small values of the arguments n, fi and v. As these arguments increase, the 
resulting expression becomes much more unwieldy. The expression simplifier 
makes use of the power of Mathematica's symbolic manipulation functions 
to obtain the most simplified expression for I(u,n, ^,1/). This task is com- 
plicated by the branch-cuts of the complex logarithm function. Due to this 
and the complexity of the starting expression, direct simplification does not 
seem to yield the most simplified expression. Rather we use expression tree 
parsing and manipulation that consists of the following steps: (i) simplify 
each Log expression in terms of the three elementary Logs that constitute 
the entire Log dependence: Log[i2M + 1], Log[i2-u — 1] and Log[-u]. This step 
is not as trivial as it seems due to the need to handle the complex arguments 
properly, (ii) Check the coefficients of the Log[i2M -|- 1] and Log[i2M — 1] are 
the same and real (as I{u, n, /i, u) is known to be real) and combine them 
into Log[4M^ + 1] + Log[— 1], and (iii) finally check that the whole expression 
is real. The Mathematica functions that are most important for such parsing 
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and simplification are the Numerator^ Denominator^ Extract, Together and 
Head. For the actual implementation, refer to the submitted Mathematica 
code. To demonstrate the power of this technique, taking —1, 1, 1) as 



given by Eq.(16) as input, the expression simplifier produces the following 
expression 

/(u, -1,1,1) = |^4m2(-1 + 6^2) - 32^3 ArcTan[2M] 

+ {I + l2u^)l.og[l + Au^]^ . (17) 
As a final test, we note that the expression simplifier indeed recovers the 



simplified results for hi{u) and h2{u) (starting from Eqs.(14) and (15)) which 



are also one obtained analytically in appendix 7.1 In the next section we 



discuss the final part of the submitted Mathematica code. 

3.3. Expression checker 

From the implementation point of view, this is the simplest part of the 
code. The expression checker tests the accuracy of the simphfied analytical 
expressions for I{u,n,fi, v) against a high precision numerical calculation of 
the same quantity. Specifically, we make used of the NIntegrate function in 
Mathematica with a precision of 10 decimal digits and compare the output 
of the symbolic with the numerical computation for a range of values of u. 



Hence, for the sample analytical expression given in Eq. (17), the numerical 
counterpart is produced by 

templist = N Integrate [Exp [-X / tempu] x" (-1) SphericalBesselJ[l, X ] 

SphericalBesselJ[l, x], {x,0,1000*tempu\ , PrecisionGoal-^ 10] , 

where templist is the value of the numerical integration for some u = tempu. 
We discuss the result in section 14.31 



4. More on the Mathematica package 

There are two files submitted to the electronic library. The first file is the 
actual Mathematica package SymbBesselJInteg.m that contains the imple- 
mentation of symbolic integration algorithm and related accuracy checking 
function. In the second file, we have implemented an extensive test of the 
algorithm. 
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4.I. Program structure 

There are three main functions made accessible in the package. These 

are 

(i) ParseCalcIInt[u, n, fi, v] which takes four arguments that correspond 
to the four variables in Eq. (|2]). The function calculates and returns 
I{u, n, II, v) without simplification, in essence implementing the expres- 



sion parser discussed in section 3.1 



(ii) Simplifier[u, lint] whose second argument should be the exact expres- 
sion returned by the ParseCalcInt function and while the first one is the 
basic variable of the problem u. The function essentially implements 



the expression simplifier discussed in section |3.2[ Its return value is the 
most simplified expression corresponding to the input expression, 
(iii) CheckAnalyNumer[u, n, fi, v, uList, lint] is the function that tests 
the accuracy of the symbolic integration algorithm against a numerical 
result of the same quantity. It takes six arguments. While the mean- 
ing of the first four arguments should be self-explanatory, the fifth one, 
uList, refers to a list containing u values for which the test is going 
to be carried out. The last argument, lint, refers to either the output 
of the ParseCalcInt function or that of Simplifier, depending on which 
code segment is being tested. In short, this function implements the 



expression checker discussed in section |3.3[ Finally, the function re- 
turns a list containing both the real and imaginarjjj parts of both the 



analytical expression and the corresponding numerical result. 

In the following, we illustrate how one can make use of the three functions 
in the package to calculate Eq.(|2]) and to verify the accuracy. 

4-2. Illustrative example 

In order to demonstrate how to use the package, we calculate I{u, —3, 3, 3). 
First, one needs to load the package into the Mathematica session. Assuming 
the package is in the "temp" directory, this is accomplished by 

In[l]: = SetDirectory["/temp"]; 
In[2]: = Needs["SymbBesselJInteg"']; 



^Thc imaginary part is expected to be zero. 
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To get a list and description of the functions made available in the package, 
one can use the command 



In [3]: = ?SynibBesselJInteg 



In the next step, one should clear the basic symbol which we call u (note 
that any variable name is valid except that one has to use it consistently) 
and call the ParseCalcInt function as 



In[4] 
In [5] 
In[6] 



— Clear[u, n, /i , u , uList]; 
= n = -3; n = 3; v = 3; 

— lint — ParseCalcInt [u ,n, /i, v]; 



At this point, the variable lint has the complicated analytical expression for 
I{u, —3, 3, 3). The simplification proceeds by invoking 

lii[7]:=IintSimpUfied = Simplifier[u ,Iint] 

Thus lintSimplified contains the simplified expression. For the case at hand, 
viz, I{u, —3, 3, 3), it reads 



Out [7]: 



Au^ (-15 - 240^2 - 1556m^ + 4272m^ + 672m^) 



967680iiio 

-1536 J (15 + 4:u^) ArcTan[2K] + 3(5 + 90^^ + 672u'^ + 3360u^) 



X Log [l + 4^2] 



(18) 



Finally, to check the accuracy, one starts by definining a list of values of u 
as, for example, 

In[8]: = uList = {0.5, 1.0, 2.0, 4.O, 8.0}; 

and then issue the command 

In[9]: — CheckAnalyNumer[u, n, /i, v, uList, lint] , 

which returns a list of values containing the results of the symbolic algorithm 
and numerical calculations. This is discussed in the next section. 
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4-3. Validity of the results 

The extensive test file is one of tlie two files submitted to the electronic 
library. In this test which is carried out over the practically important part of 
the problem domain, we did not encounter a single case where the algorithm 
and its Mathematica implementation fail to hold. This extensive test was 
necessitated by both the present work and the application mentioned in Ref. 
|15j . In addition, as n, fi and u increase, the instability of the naive numer- 
ical computation becomes too prohibitive to use. This makes the symbolic 
integrator discussed in this paper and other numerically stable packages dis- 
cussed in Refs. [101 IH] indispensable. Finally, we conclude this section with 
Table 1, that shows the comparison of the analytical and numerical values 
for I{u, —3, 3, 3) for a selected set of u values. In both cases, the imaginary 
part is zero, hence it is not shown. As can be read from the table, the two 
values match exactly at least up to the reported number of decimal digits. 





Tabic 1: I{u, 


-3,3,3) 


u 


Analytical 


Numerical 


0.5 


0.0000216 


0.0000216 


1.0 


0.0001523 


0.0001523 


2.0 


0.0005597 


0.0005597 


4.0 


0.0011940 


0.0011940 


8.0 


0.0017990 


0.0017990 



5. Comments and conclusions 

A symbolic algorithm for the analytical integration of a product of two 
spherical bessel functions with an additional exponential and polynomial fac- 
tor was discussed and implemented in a Mathematica package. The package 
complements existing numerical packages [TUl ITT] . It will especially be useful 
in situations where one has to calculate a huge number of these integrals 
analytically, as was the case in the application of the density matrix ex- 
pansion [151 113 HI] to the HF energy obatined from chiral EFT two- and 
three-nucleon interactions at N^LO. Possible directions of research could be 
to expand the package to treat general bessel functions, to allow for the ad- 
dition of a scaling parameter in the arguments of the bessel functions and to 
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handle more than two bessel functions as is done in the numerical approach 
ofRef. HI]. 



6. Acknowledgements 

We thank R. J. Furnstahl for useful discussions. This work was supported 
in part by the National Science Foundation under Grant No. PHY-0758125 
and the UNEDF SciDAC Collaboration (DOE Grant DE-FC02-07ER41457). 



7. Appendix 

7.1. Analytic evaluation of hiiu) and h2{u) 



Let us calculate the limits given in Eqs.(19) and (20) 

y 



hi{u) 
h2{u) 



lim 
lim 



Ei{ -i2y- - ] + Ei{ i2y 

u I V u 



2Ei 



u 



Ei \ -i2y 



y 



u 



Eili2y- - 



, (19) 
(20) 



Starting with the definitions of the exponential integral function given in Eq. 
(|6]) and given its analyticity over the complex plan^ performing a change 
of variable followed by a few simplification steps yield 



hi{u) 
h2{u) 



lim 

y^oo 

lim 

y^oo 



-y/u 



i2y 



lo t + y/u 
/ ——dt 

Ji2y t + y/U 



p-- , i'iy p-t 

dt + e~'^/" / -dt 



t + y/u 



(21) 
(22) 



The next several steps involve change of variable of integration, applied uni- 
formly to both hi{u) and h2{u) through rotation by 7r/2 and vr consecutively. 
Hence 



hi{u) 
h2{u) 



lim 

i/->oo 

lim 



-y/u 



2y 



2t cos t 



+ 



2y sin t 



-I e 



-y/u 



\t^ + {y/n)^ u{t^ + {y/uY) 

2s/ pit 



dt 



2?/ 



-it + y/u 



dt 



(23) 
(24) 



^It has a branch cut along ] — oo, 0]. 
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Hence, taking the limit of hi{u) directly and performing the remaining inte- 
gration of h2{u) (the real and imaginary parts separately), followed by taking 
the limit, we obtain 
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